Frustration of the interlayer coupling by mobile holes in La 2 - x Sr x Cu04 (x < 0.02) 
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We have studied the interlayer coupling in the antiferromagnetic (AF) phase of Sr and Zn doped 
La2CuC>4 by analyzing the spin flip transition in the magnetization curves. We find that the inter- 
layer coupling strongly depends on the mobility of the hole charge carriers. Samples with the same 
hole content as well as the same Neel temperature but a different hole mobility, which we adjusted 
by Zn co-doping, can have a very different interlayer coupling. Our results suggest that only mobile 
holes can cause a strong frustration of the interlayer coupling. 

PACS numbers: 74.25.Ha, 74.72.Dn, 75.30.Hx, 75.40.Cx 
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In La2- x Sr x Cu04 the 3D AF order is destroyed by a 
remarkably small amount of 2% holes (x = 0.02) Q. In 
this compound the S=l/2 Cu spins in the Cu02 planes 
form a 2D Heisenberg AF with a strong Cu-Cu in-plane 
exchange of J = 135 meV, and a very weak interlayer 
coupling J± ~ 10~ 5 J. Though the interlayer coupling 
is indispensable to establish the 3D AF order j2|, al- 
most nothing is known about its dependence on doping 
with holes and/or static spin vacancies. For La^CuC^ 
Thio et al. have demonstrated that Jj_ can be determined 
from a spin flip (SF) for H \\ c due to the Dzyaloshinsky- 
Moriya (DM) Cu-spin canting ||. 

In recent years, great efforts have been made to map 
each stage of the suppression of the 3D order which is 
a precondition for the occurrence of superconductivity 
(SC) @, § |, § . The Neel temperature T N decreases from 
325 K for x = to about 80 K for x = 0.019 and then 
drops to zero within a Sr range of Ax ~ 0.001 around x = 
0.02 H, ||. At the same time for x > 0.008 the so called 
spin freezing regime (SFR) evolves at T < 30 K ||, ||. 
At x — 0.02 the SFR crosses over to the cluster spin glass 
phase which reaches into the SC phase that appears at 
x = 0.06 §. In the AF phase for 30K < T < T N ^SR 
and NQR probe a rapid decrease of the internal field 
Hi n t with increasing x. However, in the SFR iTnt again 
rises and independent of x reaches almost the value in 
La2Cu04. Recent neutron diffraction data prove the SFR 
to have a incommensurate AF spin modulation, which 
coexists with the commensurate AF order . 

Various models were suggested to describe the sup- 
pression of the AF order §, |, |, |(], [ll). The frus- 
tration model 0] assumes that individual localized holes 
cause a frustration of the in-plane and the interlayer cou- 
pling. In the finite size scaling model holes segregate 
into domain walls which limit the 2D correlation length 
£,2D oc vT/v". The SFR was associated with the breakup 
of the domain walls when holes localize at low T In 
contrast, recent neutron diffraction data support the ex- 
istence of magnetic domain walls in particular for the 



SFR §. 

To elucidate crucial parameters for the suppression of 
the 3D AF order in La2- x Sr x Cu04, we have studied the 
interlayer coupling in Sr and/or Zn doped samples. Zn 
doping is used to reduce the mobility of the holes as well 
as to introduce spin vacancies. From our analysis of the 
SF transition we find that only mobile holes cause a dras- 
tic suppression of the interlayer coupling. Co-doping with 
Zn recovers an interlayer coupling almost as strong as in 
pure La2CuC>4 due to the localization of holes. Surpris- 
ingly, samples with similar TV can have a completely 
different interlayer coupling, depending on the hole mo- 
bility and the number of spin vacancies. 

The DC magnetization M(H) of five polycrystals 
La2_ x Sr x Cui_ z Zn z 04 (Tab. |jj) was measured using a 
VSM {T max = 290 K, H max = 14 Tesla). Samples were 
annealed in vacuum (l/2h, 800°C), their preparation was 
described in Ref. [| 

In the LTO phase of La2Cu04 the CuOg-octahcdra 
are tilted. As a consequence, nearest neighbor Cu spins 
are slightly canted against each other due to a DM su- 
perexchange term S. Below TV the canted moments 
Mum II c of adjacent planes are AF ordered, but can 
be flipped for H || c. Because of this moments of 
the order of ~ 10 _3 //s/Cu the M(H) curves are non- 
linear in the AF phase as well as in the paramagnetic 
(PM) phase || This can be seen in Fig. |l| where 

we compare M(H) of La 2 Cu0 4 (T N = 312 K) and 
Lai.983Sro.oi7Cu04 (TV = 132 K) as a function of T/T/v. 
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312 


4.5(5) 


2.7(3) 


2.9(5) 


0.011 





222 


3.6(3) 


2.1(3) 


1.7(5) 


0.017 





132 


2.4(2) 


1.2(2) 


0.7(3) 


0.017 


0.10 


134 


3.9(4) 


3.0(3) 


2.7(5) 





0.15 


166 


3.6(4) 


3.4(3) 


2.7(5) 



TABLE I: Studied La2- x Sr x Cui_ z Zn z 04 samples (see text). 
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FIG. 1: Magnetization curves M(H) of La2- x Sr x Cu04 with 
x = and x = 0.017 at reduced temperatures T/Tjv = 1.05, 
0.9, and 0.2 (see text). Inset: M(14T) vs T/T N . For x = 
data only for T/Tjv < 1 as T N > T max (VSM). 
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FIG. 2: Spin flip term Msf = M - Mb -\oH in La 2 Cu0 4 at 

30 K. ( )/ ( ) fit according to Eq. H without/with gaussian 

distribution of Hsf- Left inset: Spin flip takes place when 
HcosO \\ c exceeds Hsf- Only MpfjCosO || H contributes to 
Msf- Right inset: Derivative cLMsf / dH . 



In the PM phase (T/Tjv > 1) we find that the non-linear 
contribution can be described by a Brillouin term Mb- 
In the AF phase (T/T/v < 1) the non- linear contribution 
is always a combination of Mb and a term Msf which 
arises from the SF of the DM moments: 

(xoH + M B (T > Tjv) m 

\ xoH + M b + Msf {T<T n ) 1 1 

where xoH accounts for all linear terms. As we will see 
below, Msf is a function of the long range AF ordered 
fraction of the DM moment M^j(T) cx AF order pa- 
rameter, and Mb is a function of the non-AF-ordered 
(NO) fraction M^°, both per Cu spin. Just below T/v 
M(H) for x = and x = 0.017 is almost identical up to 
14 Tesla and Msf is very small (grey shaded) while Mb 
is still large. In contrast, at low T the non-linear part of 
M(H) strongly depends on x. For x = it is dominated 
by Msf whereas Mb — 0. For x = 0.017 Msf is much 
smaller and of about the same order as Mb- Clearly vis- 
ible is also a much larger Hsf for x = 0. At 14 Tesla M 
is almost x-independcnt for all T/Tjv (inset Fig. 1). 

The analysis of the SF enables to extract Hsf(T) and 
m dai( t )- From the low T limes J J- can be calculated [§: 

M^(0)H SF (0) ~S 2 J ± . (2) 

The SF takes place when H acting on Mp^ (see left 
inset Fig. ^ overcomes the interlayer coupling, i.e. for 
HcosO > Hsf \\ c all M^f { become ferromagnetically 
aligned || c. In a single crystal (H || c) the SF causes 
a step like increase of M(H) at H ~ Hsf by Msf = 
Mpj^. In a polycrystal the crystallites are oriented ran- 
domly. If Hs f is identical in all crystallites we obtain the 
following field dependence of Ms f for H > Hs f from in- 
tegration over all directions: 

Msf(H) = l/2Af££[l - (Hsf/H) 2 }. (3) 



Obviously in a polycrystal Msf converges to l/2M^ f 
for H — > oo. In Fig. ||we show Msf = M— Mb— XoH for 
La2Cu04 at 30 K. For large H Eq. || (dashed line) yield 
a good fit to the data. Assuming a gaussian distribution 
of Hsf in the polycrystal, integration of Eq. || over Hsf 
yields the solid line in Fig. || which perfectly fits the data. 
The extracted parameters M^ t = 2.55 x 10 -3 fXs and 
Hsf = 5.2 Tesla are in fair agreement with single crys- 
tal data in Rcf. ||. Note, that in a polycrystal Hsf is 
located slightly below the maximum of dM/dH (inset). 

For the Brillouin term Mb in the PM phase of a poly- 
crystal a similar geometrical consideration as for Eq. || 
yields the phenomenological formula: 

mt/2 

M B (H) = Mqm I tanh(kHsin<t)) sin 2 (j> d(f>. (4) 
Jo 

where k = M%$N/(k B T + Jj_N 2 S 2 ) is a phenomeno- 
logical expression with N = (£215 / a ) 2 the number of 2D 
correlated Cu spins, which provides a rough estimate for 
t;2D ■ <j> is the angle between H and the octahedra tilt axis 
which is normal to the DM plane (cf. inset Fig. |J). In the 
PM phase M$& equals the full DM moment M DM and 
for H — ► co Mb converges to tt/4Mdm- With k used as 
a simple fit parameter Eq. ^ describes also the remain- 
ing contribution Mb of the non-AF-ordered fraction of 
Mbm in the AF phase. Our analysis has shown that 
M DM ^ M£%j + M%£. A neglect of M B would lead to 
inaccurate values for M^^i Hsf, and J±, in particular 
when M B > M S f- 

In Fig. U we show the T dependence of M^^ and Hsf 
of the five samples studied (cf. Tab. Q). In pure La2Cu04 
M^M increases monotonous with decreasing T, and ex- 
trapolates to A^m(°) = 2 - 7 x 10_3 Mb/Cu for T -> K. 
Pure Sr doping causes a drastic reduction of M^^. In 
contrast, in La2Cuo.85Zno.15O4, though Tjv is strongly 
reduced, at low temperatures M^m (per Cu atom) be- 
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FIG. 3: Spin flip parameters M^m ( a ) an d Hsf (b-e) vs 
T in La 2 _xSr x Cui_ z Zn z C>4. (b-d) H SF for dH/dt > and 
dH/dt < as well as mean value (x). (e) Hsf mean values. 

comes even larger than in La2Cu04. Interestingly, a sim- 
ilar behavior is observed for the Sr/Zn co-doped sample, 
where the twofold role of Zn is to create spin vacancies 
and to reduce the mobility of the holes. 

In La2Cu04 Hsf increases with decreasing T and be- 
low 150 K becomes hysteretic, with its mean value satu- 
rating at 4.5 Tesla [Fig. ||(b),(e)] . Pure Sr doping strongly 
reduces Hsf as well as the hysteretic T range. In partic- 
ular for 1.7% Sr only for T < 10 K a clear field hysteresis 
is observed. Both Zn doped samples show a relatively 
large field hysteresis [sec (c),(d)]. Their maximum mean 
value for Hsf is smaller than in La2Cu04, but in view 
of their relatively low Tjv, Hsf is large if compared to 
Lai.983Sr .oi7Cu04. 

As one can see in in Fig. |^(a), in La2CuC>4 and for 
pure Sr doping decreases below 30 K. Though 

this looks similar for all x, there is a qualitative differ- 
ence. In La2CuC>4 the hysteresis of the SF increases up 
to 2.2 Tesla at low T and M(H) becomes strongly dis- 
torted which causes deviations from the fit and makes 
it difficult to extract M^y. In contrast, in particular 
in the 1.7% Sr doped sample drops between 30 K 

and 15 K where M(H) is reversible. Therefore, we think 
that here the drop of M^jfj is connected to the SFR. 
As Mjjfj oc AF order parameter, its drop in the SFR 
signals a degradation of the long range order, consistent 
with the neutron diffraction data in Ref. 0. Obviously, 
this effect is absent in both Zn doped samples. For the 
Sr/Zn doped sample this evidences for a suppression of 
the SFR in La 2 - x Sr x Cu0 4 by Zn. 

Following Eq. we show in Fig. [| the T-dependence 
of J* = M^jHsf/S 2 which we call the effective in- 
terlayer coupling, where Hsf are the mean values in 
Fig. ||(c). J* accounts for the effects of doping (x,z) 
and temperature which counter the interlayer superex- 
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FIG. 4: Interlayer coupling in La2- x Sr x Cui_ z Zn z 04. Inset: 
Inverse correlation length vs T for both samples with x = 
0.017. Solid lines are guides to the eye. 



change and only for T -> J£(p) = J± (sec Tab. |). In 
La2Cu04 Jl_ increases with decreasing T and for T — > 
we find JTl(0) = 2.9 ^eV which is in good agreement 
with 2.6 ficV for a single crystal with Tn — 240 K 
in Ref. As a function of Sr doping drastically 
decreases, and for x = 0.017 J7^(0) amounts to only 
25% of the value in La2Cu04. If we compare this to 
La1.983Sro.017Cuo.9Zno.1O4 which has exactly the same 
Tjv, the huge difference is apparent. Here, J± rises 
steeply and reaches almost the same value as in La2Cu04. 
La2Cuo.85Zno.15O4 behaves similar. With the hindsight 
that i7J(0) in LaaCuC^ and the Zn doped samples is 
about the same, it is evident from Eq. || that the larger 
DM moments in the Zn doped samples have to result in 
correspondingly smaller critical fields Hsf [Fig. ||(a),(e)]. 
To a certain extend the 15% (z = 0.1) and 25% {z = 0.15) 
larger DM moments can be explained by an enhanced 
octahedra tilt angle $ [Q . Local strain around the non- 
Jahn- Teller-active Zn sites might amplify this effect. 

From Fig. |] it is obvious that the mechanism for the 
suppression of the 3D AF order is completely different for 
pure hole doping (Sr) on the one hand and non-magnetic 
impurity doping (Zn) as well as mixed Sr/Zn doping on 
the other hand. In La2- x Sr x Cu04 the hole mobility in- 
creases rapidly with increasing x, and just a small con- 
centration of mobile holes reduces both, Tn and J~*. Co- 
doping with Zn reduces the hole mobility || and causes 
a drastic increase of , even if Tn is low. In particular 
for the two samples doped by 1.7% Sr the influence of the 
hole mobility is apparent. Our data clearly demonstrate 
that at fixed Sr content x < 0.02 the interlayer coupling 
drastically increases if the hole mobility is reduced. 

Localized holes are supposed to suppress Tn much 
stronger than Zn because a hole located on an O-site fer- 
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romagnetically frustrates the AF Cu-O-Cu exchange 
In contrast, a static spin vacancy does not perturb the 
AF order of the surrounding Cu spins E3]. If we com- 
pare now the 10% Zn and 1.7% Sr co-doped sample with 
a ~18% Zn doped sample, which would have the same Tv 
of ~T35 K, it is evident that 1.7% localized holes suppress 
Tv as effective as 8% Zn. This means that also for local- 
ized holes the effect of in-plane ferromagnetic frustrations 
on Tjv is strong, while their effect on the interlayer cou- 
pling is very weak (Fig. ^), which is in contrast to the 
frustration model 0. Obviously the suppression of the 
interlayer coupling in La2_ x Sr x Cu04 is closely connected 
to a high hole mobility (cf. Ref. |(] for p values). 

Lets focus on the correlation length of the two samples 
with 1.7% Sr in the inset of Fig. [|. According to our 
analysis they have about the same ^2D at Tv- However, 
the sample with localized holes shows a much faster in- 
creases of £-2D with decreasing T than the sample with 
mobile holes. In view of these results for T > T/v and the 
fact that the AF ordered moment at low T stays much 
smaller in the purely Sr doped sample, we think that 
mobile holes hinder ^2D to develop and frustrate the in- 
terlayer coupling in a way that cannot be explained with 
isolated localized holes. 

Most of our observations presented in this paper can 
be naturally explained assuming dynamic magnetic an- 
tiphase boundaries for T > 30 K and static for T < 30 K. 
Evidence for static antiphase boundaries in the SFR was 
recently found by Matsuda et al. R]. Below 30 K an 
incommensurate AF order was detected which coexists 
with the commensurate AF order. Dynamic antiphase 
domains were suggested to explain the drastic reduc- 
tion of the AF order parameter with increasing x for 
30 K < T < T/v jn|, but so far there is no direct ev- 
idence for that. Our discussion does not depend on a 
particular domain form or on the formation of charge 
stripes. We think that what is essential is the presence 
of mobile holes, as only these holes can affect many Cu-O- 
Cu-bonds and are able to excite at least antiphase bound- 
ary segments, i.e. strings of ferromagnetically frustrated 
Cu-O-Cu-bonds. In such an environment a frustration 
of the interlayer coupling is inevitable, as two antiphased 
domains in one Cu02-planc cannot AF couple simultane- 
ously to the adjacent plane, assuming that in particular 
in the case of dynamic antiphase boundaries 3D correla- 
tions are absent. 

The drastic reduction of M^f { in La2- x Sr x Cu04 with 
increasing x for 30K < T < Tv as well as its further 
reduction in the SFR can be explained in this framework 
too. Just those regions of the Cu02 planes with AF in- 
terlayer coupling contribute to the spin flip (Msf)i while 
regions with ferromagnetically frustrated interlayer cou- 
pling can be aligned continuously with increasing field 
(Mb)- Both contributions must result from correlated 
regions with relatively large ^2D as it is possible to align 
the full DM moment at any temperature (cf. Fig. 1 for 



high H up to H = 14 Tcsla). Interestingly, no signif- 
icant change of the mean value H$f is observed upon 
the transition into the SFR, whereas the decrease in 
M$f and a corresponding increase in Mb was very clear. 
Obviously, just the fraction of regions with AF inter- 
layer coupling decreases and not their interlayer coupling. 
This suggests that in the SFR coupled and frustrated re- 
gions are spatially separated, which is consistent with 
the observation of two /j,SR frequencies in the SFR of 
the Lai.g83Sro.oi7Cu04 sample |l5| and recent neutron 
diffraction data in Ref. |?]. The reduction of in the 

SFR could be the consequence of a more destructive or 
a more even distribution of antiphase boundaries when 
they become statically ordered. We think that on either 
side of an in-plane antiphase boundary the DM moments 
point into opposite directions. At first glance this looks 
like DM moments cancel out, but it is important to con- 
sider both, their arrangement in the planes as well as be- 
tween the planes. In-plane phase boundaries cause jumps 
in the AF in-plane phase but do not necessarily destabi- 
lize in-plane correlations, whereas they certainly result in 
a frustration of the interlayer coupling. 

On the basis of our experimental data we come to the 
conclusion that the drastic suppression of the interlayer 
coupling and hence of Tv in La2- x Sr x Cu04 is a real mag- 
netic decoupling of the Cu02 planes due to the ferromag- 
netic frustration of the interlayer coupling by in-planc 
magnetic antiphase boundaries. We have shown that an 
important condition for this effect is a high hole mobility, 
and that the ferromagnetic frustration effect of localized 
holes on the interlayer coupling is small. 
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